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(57) A method for removing "grid line artifacts" from 
x-ray .mages without changing the diagnostic quality of 
the x-ray image is presented The method utilizes the 
Founor spectrum (216) of the image to detect the grid 
■■re frequences and employs spectral domain filtering 
to remove the gr,d line spectral components. The diag- 
nostic information is preserved by modifying (214 212) 
the g-.d line spectral components to be indistinguishable 
from local variations in image intensity values and edqe 
dons ly of the x-ray image 
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D scription 

[0001] This invention relates to a Fourier spectrum 
method for removing "grid line artifacts" from x-ray im- 
ages, and more particularly to a Fourier spectrum meth- 
od for removing "grid line artifacts" without changing the 
diagnostic quality in x-ray images. 
[0002] In a x-ray radiography imaging system, an x- 
ray source projects a cone-shaped pattern of beams. 
This cone beam passes through the object being im- 
aged, such as a medical patient, and impinge upon a 
two-dimensional array of radiation detectors. The signal 
generated from the measurement of the intensity of the 
transmitted radiation is dependent upon the attenuation 
of the x-ray beam by the object. Each detector produces 
a separate electrical signal that is a measurement of in- 
cident beam attenuation. 

[0003] A metallic anti-scatter grid used in the x-ray ra- 
diography imaging system is typically placed against the 
detector array to allow the x-rays that trend along a sub- 
stantially perpendicular path to the a respective detector 
to strike the respective detector, and x-rays that do not 
trend along a substantially perpendicular path to the de- 
tector are blocked by the anti-scatter grid, as is illustrat- 
ed in Figure 1 As such, the anti-scatter grid enhances 
image diagnostic quality by preventing undesired x-rays 
from striking a detector. A disadvantage of using the anti 
scatter grid is that it may cause "grid line artifacts" to 
appear in the x-ray image. The grid line artifact appears 
in x-ray images as intensity modulation of the image in 
lines parallel to the anti-scatter grid. The grid line artifact 
occurs when the grid lines run perpendicular to the scan 
lines on the display device. The grid line artifact is very 
sensitive to display image magnification and can be 
made worse or be made to disappear by changing the 
image magnification. It would be desirable to remove the 
"grid line artifacts" from the x-ray image without chang- 
ing the diagnostic quality in the x-ray image, irrespective 
of image magnification. 

[0004] Another cause of "grid line artifacts" on the x- 
ray image include errors generated by the detector array 
signal measurement electronic circuits. For example, 
when there is a difference in gain between two respec- 
tive signal measurement circuits, "grid line artifacts" 
may appear in the x-ray image. 

[0005] Yet another source of "grid line artifacts" may 
be caused by the repositioning of the anti-scatter grid at 
non-standard positions during successive scans called 
"over-sampling. " Additionally, "grid line artifacts" may be 
caused by variations in the x-ray dose during over-sam- 
pling. It would be desirable to remove grid lines from the 
x-ray image caused by over-sampling without changing 
diagnostic quality. 

[0006] A method for removing "grid line artifacts" from 
x-ray images in an x-ray radiography imaging system is 
presented. The method utilizes the Fourier spectrum of 
the image to identify grid line frequencies and employs 
spectral domain filtering to remove the grid line spectral 


components. The diagnostic information is preserved by 
modifying the grid lin spectral components so as to be 
indistinguishable from local variations in image intensity 
values, and edge density of the x-ray image. Grid line 

s spectral components are removed by the method con- 
sisting of the following steps: first, replacing the edgy 
regions with non edgy regions to generate a modified x- 
ray image; next, replacing the high intensity regions with 
low intensity regions within the modified x-ray image; 

io next, converting the modified x-ray image to the fre- 
quency domain; then, eliminating "grid line artifacts" 
from the modified x-ray image; and finally, converting 
the modified x-ray image to a human readable format. 
[0007] The features of the invention may best be un- 

75 derstood by reference to the following description in 
conjunction with the accompanying drawings in which 
like characters represent like parts throughout the draw- 
ings, and in which: 

[0008] Figure 1 is an illustration of x-rays bombarding 
20 a prior art anti-scatter grid and the detector array of a x- 
ray radiography system. 

[0009] Figure 2 is a method flowchart of the method 
for removing "grid line artifacts" of the present invention. 
[0010] Figure 3 is an illustration providing further de- 
2S tail of the method for removing "grid lino artifacts" of the 
present invention. 

[0011] Figure 4 is a continuation of the method for re- 
moving "grid line artifacts" of Figure 3. 
[0012] Figure 5 is a continuation of the method for re- 
30 moving "grid line artifacts" of Figure 3 

[0013] Figure 6 is a further illustration of block 2B0 il- 
lustrated in Figure 4. 

[0014] Figure 7 is an illustration of the spectral com- 
ponents of a Fourier transform of the typical x-ray image 

35 with "grid line artifacts". 

[0015] Figure B is an illustration of spectral compo- 
nents of the Fourier transform illustrated in Figure 7 hav- 
ing "grid line artifacts" modified by the method of the 
present invention. 

40 [0016] In a x-ray radiography imaging system, an x- 
ray source (not illustrated) projects a cone-shaped pat- 
tern of beams which pass through an object 124 being 
imaged and impinge upon an array of radiation detec- 
tors 130, as illustrated in Figure 1. Each detector 130 

45 produces a separate electrical signal that is a measure- 
ment of incident beam attenuation. The attenuation 
measurements from all detectors 1 30 are acquired sep- 
arately to produce an x-ray image. A metallic anti-scatter 
grid 114 is typically placed against a detector base 112 

50 to channel the x-rays so that only substantially perpen- 
dicular x-rays (e.g. rays 118 and 122) to a detector strike 
the detector 1 30 and x-rays that are not substantially 
perpendicular (e.g. x-rays 116 and 120) arc blocked by 
anti-scatter grid 114. Substantially perpendicular x-rays 

55 are those x-rays that strike detector 1 30 and do not strik- 
ing anti-scatter grid 114. 

[0017] In this Specification "grid line artifacts" are de- 
fined as extraneous visible data in Ihe x-ray image gen- 
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erated in association with the use of anti-scatter grid 
114, detector 130 electronics, and any other source 
which causes extraneous data to be seen in the x-ray 
image which result in high frequency spectral compo- 
nents having a spectral magnitude significantly greater 
than high frequency object data spectral components. 
Object data spectral components comprise those fre- 
quency components that result from the x-ray image and 
not from artifacts. 

[0018] The method of removing "grid line artifacts" is 
illustrated by the method flowchart shown in Figure 2. 
This method includes: a step 210, subdividing the x-ray 
image into a plurality of windows; a step 212, replacing 
regions with a substantial number of edges within each 
window with regions with less number of edges in the 
respective window; a step 214, replacing high intensity 
regions within each window with low intensity regions; 
a step 216, converting each window into the Fourier do- 
main and removing "grid line artifacts"; and a step 216, 
conducting an inverse-Fourier transform on each win- 
dow so that the modified x-ray image may be presented 
in a human readable format. Regions with a substantial 
number of edges hereinafter wilt be identified as "edgy 
regions" in this Specification. Edgy regions are further 
discussed below. 

[0019] The x-ray image is divided into windows to fa- 
cilitate computer based evaluation, as identified by the 
method of step 210 in Figure 2. The window size is se- 
lected to be large enough to generate frequency reso- 
lution at about the expected frequency of the grid line 
artifact and alternatively, small enough to insure that the 
grid line artifact has a substantially constant spatial grid 
line frequency. The size of each window selected is such 
that it has substantially constant grid frequency. A grid 
frequency is considered substantially constant if there 
is at most one local maximum frequency component 
which is three standard deviations larger than the aver- 
age value of the neighboring frequency components. 
For example, in a typical x-ray image having 1024 by 
1024 pixels, the image may be divided into sixty-four 
windows. It is also understood that a window may have 
rectangular dimensions rather than square dimensions 
as illustrated in this case. 

[0020] The "edgy region" of the x-ray image provides 
the necessary definition of object 124 (Figure 1) so that 
the human observer may easily identify important fea- 
tures of the x-ray as identified by the method of step 2 1 2 
in Figure 2. These features aid in diagnosis. For this rea- 
son it is important to preserve the detail of the "edgy 
region" of the x-ray image. The method employed in this 
invention preserves the "edgy regions" of the x-ray im- 
age by removing them before processing the image and 
replacing thorn after processing the image. "Edgy re- 
gions" are defined by the method identified below. 
[0021] High intensity regions of the x-ray image are 
also removed before processing, as identified by the 
method of step 214 This is necessary first to preserve 
the high intensity regions produced by the object and 
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second, to facilitate removal of "grid line artifacts 1 ' in win- 
dows with high intensity regions. 
[0022] Next, "grid line artifacts" are removed from the 
x-ray image, as identified by the method of step 216. 

s This is accomplished by converting the x-ray image to 
the Fourier domain where frequency components rep- 
resenting the x-ray image are generated, as illustrated 
by graph 370 in Figure 7. in Figures 7 and 8 the vertical 
axis is defined as the frequency component magnitude, 

10 expressed in Decibels. The horizontal axis is defined as 
the frequency axis, expressed in Hertz. The frequency 
component of a grid line artifact 380 is distinguished 
from the frequency components of the image (e.g. com- 
ponents 372 - 378, and 382 - 386) in one respect in that 

is spectral component 380 is usually the highest spectral 
component of the high frequency components. As such, 
by eliminating low frequency components, (i.e. those 
components identified by range 371), the grid line arti- 
fact is usually the highest spectral component remain - 

20 ing, as is further discussed below. Grid line spectral 
component 380 is adjusted so that component 380 is no 
greater in magnitude than the average adjacent magni- 
tude. The adjacent magnitude is defined by range 388 
as illustrated in Figure 7. 

2S [0023] The x-ray imago is then restored to a human 
readable format, as identified by the method of step 218 
of Figure 2. This is accomplished by conducting an in- 
verse Fourier transform on the frequency components. 
It may also be necessary to adjust the mean value of 

30 the image, as is discussed below 

[0024] Method 200 for removing "grid line artifacts" as 
shown in Figure 2, is presented in greater detail in meth- 
od 250, as illustrated in Figure 3 through Figure 5. The 
steps identified in Figure 3 identify the method of f gen- 

35 erating a gradient image of the x-ray image and an in- 
tensity quantized image from the x-ray image. The gra- 
dient image is generated by first convolving the image 
with a gradient operator, as identified by step 254. For 
digital images, such as those in this invention, gradient 

40 operators represent finite difference approximations of 
either the orthogonal gradient or the directional gradient. 
These masks are defined, for example, in the reference 
entitled, Fundamentals of Digital Image Processing , Anil 
K. Jain, Prentice Hall, 1989, chapter 9, pages 347 

45 through 350. 

Next, a histogram of the gradient image is computed, as 
identified by step 260. From the histogram a threshold 
"T" is chosen such that a range of the pixels with the 
largest gradient are defined as edges, as identified by 

50 step 262. For example, threshold value "T" may be se- 
lected from a range to identify about 1 0% to about 20 % 
of the pixels with the largest gradient. A gradient image 
is therefore generated using the threshold T" to define 
edges, as identified by step 266 

55 [0025] An intensity quantized image is generated by 
first computing an intensity histogram of the image, as 
identified by step 270 in Figure 3 From this histogram 
a threshold "l" is chosen so that a fixed percentage of 
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the pixels with the larg st intensity ar defined as high 
intensity regions, as identified by step 258. Using 
threshold T, an intensity quantized image is generated, 
as identified by step 264. Regions with a low intensity 
are associated with a zero and regions with a high in- 
tensity are associated with a one. The ones and zeros 
are grouped together to define a matrix identified as the 
intensity quantized image. Intensity threshold V may be 
selected, for example, such that at most about 25% of 
the pixels with the largest intensity are defined as "high 
intensity regions. * 

[0026] The gradient image, intensity quantized image 
and x-ray image are all then divided into a plurality of 
number of respective windows as identified by step 210. 
[0027] Next, the mean value of each window is com- 
puted as identified by step 274, in Figure 4. It is neces- 
sary to determine the mean value ol the pixels of each 
window because mathematical manipulation as defined 
in subsequent operations within this Specification may 
change Ihe mean value ol the image within a respective 
window, resulting in undesirable artifacts. Thus, it may 
be necessary to adjust the mean value of each respec- 
tive window so that the mean value of the respective 
window remains unchanged, as is discussed below. 
[0028] The next stop in method 250 is to remove high 
frequency edges and high frequency intensity regions 
as identified by the method identified in step 280, in Fig- 
ure 4. In order to remove "grid line artifacts" from the x- 
ray image it is first necessary to remove high frequency 
edges and high frequency intensity regions from each 
window, as discussed above. Block 260 is further sub- 
divided into several sub-steps, as is identified by the 
method in Figure 6. 

[0029] First, a window comprising a respective inten- 
sity window, gradient window, and x-ray window is sub- 
divided into a plurality of sub-windows so that most sub- 
windows are substantially non-edgy, as identified by 
sub-step 282 ; in Figure 6. For example, in an x-ray im- 
age having 1024 by 1024 pixels wherein sixty-four win- 
dows are selected, each window may be divided into 
sixty-four sub-windows. 

[0030] Next, edgy sub-windows are replaced with the 
nearest non-edgy sub-windows based on the results 
from the test illustrated by the method in sub-step 284. 
Where twenty or more percent of the gradient pixels of 
a sub-window are larger than threshold T" the sub-win- 
dow is replaced by the nearest non-edgy sub-window, 
as identified by sub-step 292. Non-edgy sub-windows 
are defined as those sub-windows wherein less than 
20% of the gradient pixels are larger than threshold T\ 
The nearest sub-window is determined by comparing 
the distance between the center pixel of alternative sub- 
windows meeting the named criteria to the center pixel 
of the sub-window being replaced, then selecting the 
sub-window meeting the named criteria, and being clos- 
est in distance to the sub-window being replaced. 
[0031] Following sub-step 292. high intensity sub- 
windows are replaced by low intensity sub-windows and 


sub-windows which are non-edgy, as defined by sub- 
step 296. A high intensity sub-window is defined to be 
one in which 50% or more of the pixels have an associ- 
ated intensity quantized number of one. When a high 
5 intensity sub-window is identified, as illustrated by sub- 
step 294, the high intensity sub-window is replaced by 
the nearest low intensity sub-window, as illustrated by 
sub-step 296. A low intensity sub-window is defined as 
a sub-window where less than 50% of the pixels in the 
10 intensity quantized sub-window are associated with a 
one and sub-windows which are non-edgy, as previous- 
ly defined. This method is repeated until all sub-windows 
have been evaluated for high intensity regions and edgy 
regions, as identified by sub-step 298. 
[0032] After edgy regions and high intensity regions 
have been removed from the image, "grid line artifacts- 
are removed, as identified by steps 302 and 303, in Fig- 
ure 4. Removing "grid line artifacts" is best illustrated by 
the graphs 370 and 392 in Figures 5 and 6. Figure 5 
illustrates a typical graph 370 of the magnitudes of the 
frequency components of a Fourier transform of a re- 
spective window of the x-ray image having object data 
spectral components and grid line artifact spectral com- 
ponents. Grid line artifact spectral component 380 is dis- 
tinguished from object data spectral components of the 
x-ray image window because grid line spectral compo- 
nent 380 is typically the high frequency component with 
largest magnitude. Grid line artifact spectral component 
380 is the spectral component having a magnitude 
about three standard deviations greater than the aver- 
age magnitude of the object data spectral components 
in a range 388. Range 388 is defined as the range of 
object data spectral components having substantially 
constant magnitudes. 

[0033] Grid line spectral component 380 is removed 
by first removing the low frequency spectral compo- 
nents in a range 371 of graph 370. Range 371 is defined 
as object data spectral component range having object 
data spectral component frequencies in the bottom ten 
percent of the spectral component frequency range. 
Component 380 is then removed as replaced by a new 
spectral component 381 that has a magnitude substan- 
tially equally to the average of adjacent spectral compo- 
nents identified in range 388. This spectral component 
modification is illustrated by graph 392 in Figure 8. 
[0034] Although Figures 7 and 8 illustrate the grid line 
artifact spectral component removal method for a one- 
dimensional Fourier transform, a two-dimensional Fou- 
rier transform may be utilized to remove grid line artifact 
spectral components having "x" and w y M Cartesian coor- 
dinates. Steps in the two-dimensional method are sub- 
stantially the same as those defined above in the one- 
dimcnsional method, with the primary difference being 
that spectral components have an "x" and a "y" coordi- 
nate defining each component's location, whereas in the 
one-dimensional Fourier transform magnitudes of the 
speclral component will only have one coordinate defin- 
ing each component's location The two-dimensional 
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Fourier transform approach may be pr ferr d over the 
one-dimensional Fourier transform approach if the grid- 
line artifacts are not perpendicular to the V and V co- 
ordinates of the image. 

[0035] The step of dividing the x-ray image into win- 
dows and variations in the mean value of the image with- 
in a window may cause noticeable artifacts known as a 
"ringing artifacts" along the periphery of a number of 
windows. "Ringing artifacts" occur because image data 
utilized in the Fourier transform is not infinite. Theoreti- 
cally, a Fourier transform is a method of transforming 
data from the spatial domain to the frequency domain 
utilizing an infinite series, however, the x-ray images uti- 
lized in this invention have a finite number of windows 
and as such do not closely approximate an infinite se- 
ries. 

[0036] Variations in the mean value of the pixels within 
each window may be removed by subtracting the error 
signal generated by calculating the mean value of the 
pixels within the initial window identified in step 210, of 
Figure 2, and the mean value of the pixels of the modi- 
fied window identified in step 218, of Figure 2, and then 
adjusting each modified window by subtracting the pixel 
mean value of the initial window from the pixel mean 
value of the modified window, as identified by step 274, 
in Figure 4. The modified window is defined as the win- 
dow in which "edgy regions," "high intensity regions" 
and "grid line artifacts" have been removed by the 
above-described method 

[0037] Next, "ringing artifacts" caused in part by divid- 
ing the x-ray image into widows may be eliminated by a 
method called "window substitution." In this method 
each respective window along the perimeter of a re- 
spective modified window are replaced by a equivalent 
respective modified window to generate a super-win- 
dow An equivalent respective modified window is de- 
fined as a copy of the respective modified window. The 
resulting super-window has nine identical modified win- 
dows. Method 250 is then continued beginning at step 
302. in Figure 4. where the super-window is substituted 
for the modified window. Mathematical operations per- 
formed on the super-window causes ringing effects that 
normally would have appeared in the periphery of the 
modilied window to appear in the outer periphery of the 
equivalent respective modified windows. Each respec- 
tive equivalent modified window is discarded leaving on- 
ly the modified window having no "ringing artifacts." 
[0038] The x-ray image is restored next to a human 
readable format by performing an inverse Fourier trans- 
form converting the image from the frequency domain 
back into the original spatial domain image, as identified 
by the method of step 218. Next, all edgy regions and 
high intensity regions that were replaced, as identified 
by steps 21 2 and 21 4 of Figure 2 are added back to the 
modified image as identified by step 306, of Figure 5. 
[0039] The windows are then equalized by subtract- 
ing the pixel mean value of the modified window from 
the pixel mean value of the initial window to generate a 


difference value, and multiplying the r suiting difference 
value by each pixel in the modified image to obtain an 
adjusted modified image, as identified by step 308, in 
Figure 5. 

s [0040] Finally, the steps of method 250 described 
above are repeated for each window until all windows 
have been evaluated and processed according to the 
method described above as identified by steps 310 and 
312, in Figure 5. At the conclusion of this method "grid 
io line artifacts" have been removed from the x-ray image, 
"ringing artifacts" caused by this method have been 
eliminated, and the x-ray image has been converted to 
a human readable format for diagnosis and review, as 
illustrated by step 314, in Figure 5. 
is [0041] The present invention provides a method for 
removing •grid line artifacts" from an x-ray image. This 
method utilizes the Fourier spectrum of the image to de- 
tect grid line frequencies and employs a spectral domain 
filtering method to remove grid line spectral compo- 
se nenls. The diagnostic information is preserved by mod- 
ifying the grid line spectral components so as to be in- 
distinguishable from local variations and image intensity 
values and edge density of the x-ray image. 

25 

Claims 

1. A method of removing grid line artifacts from an x- 
ray image having edgy regions, and high intensity 
30 regions, said method comprising the steps of: 

generating a gradient image and a intensity im- 
age from said x-ray image; 
replacing said edgy regions with non edgy re- 
35 gions to generate a modified x-ray image using 

said gradient image; 

replacing said high intensity regions with low in- 
tensity regions within said modified x-ray image 
using said intensity image: 
40 converting said modified x-ray image to the 

Fourier transform domain; and 
eliminating grid line artifacts from said modified 
x-ray image. 

45 2. The method as recited in claim 1 , further comprising 
the step of eliminating ringing artifacts from said 
modified x-ray image after replacing said high in- 
tensity regions. 

50 3. The method as recited in claim 2, further comprising 
the step of restoring said modified image to a hu- 
man readable format after eliminating ringing arti- 
facts. 

55 4. The method as recited in claim 3, further comprising 
the step of dividing said x-ray image into a prede- 
termined number of respective x-ray windows, di- 
viding said gradient image into a predetermined 
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number of gradi nt windows, and dividing said in- 
tensity image into a predetermined number of inten- 
sity windows. before replacing said edgy regions 
and before replacing said high intensity regions. 

5 

5. The method as recited in claim 4, wherein the step 
of dividing said image into intensity windows, gra- 
dient windows, and x-ray windows further compris- 
es the step of dividing each respective intensity win- 
dow, gradient window, and x-ray window into a pre- io 
determined number of respective sub-windows. 


6. The method as recited in claim 5, wherein the step 
of dividing said image into intensity windows, gra- 
dient windows, and x-ray windows further compris- 
es the step of dividing each respective intensity win- 
dow, gradient window, and x-ray window into sixty- 
four respective sub-windows. 

7. The melhod as recited in claim 5, wherein the step 
of generating said gradient image comprises the 
steps of: 

convolving said x-ray image with a gradient op- 
orator so as to generate a gradient image; 
computing a histogram of said gradient image; 
and 

choosing a threshold value M T" so as to identify 
a predetermined range ot pixels above the 
threshold value "T", wherein the pixels above 
threshold value "T" are edges. 

8. The method as recited in claim 7, wherein the step 
of replacing edgy regions further comprises the step 
of selecting threshold value "T" so that an edgy sub- 
window is identified when about ten percent to 
about twenty percent of pixels within said sub-win- 
dow are above threshold value "T". 

9. The method as recited in claim 8, wherein the step 
of generating said gradient image further comprises 
the step of selecting said gradient operator to be a 
finite difference approximation of said respective or- 
thogonal gradient 

10. The method as recited in claim 9, wherein the step 
of generating said gradient image comprises the 
step of replacing each respective intensity window, 
gradient window, and x-ray window sub-windows 
with the nearest non-edgy sub-window when a re- 
spective gradient sub-window is edgy 

1 1 . The method as recited in claim 1 0, wherein the step 
of replacing edgy regions further comprises the step 
of selecting said predetermined percentage of gra- 
dient pixels to be about twenty percent. 

12. The method as recited in claim 5. wherein the step 


of replacing high int nsity regions comprises the 
steps of: 

computing an intensity histogram of said im- 
age; 

choosing a threshold value T such that a pre- 
determined percentage of pixels with the larg- 
est intensity is a high intensity region; and 
generating an intensity quantized image having 
high and low intensity regions. 

13. The method as recited in claim 12, wherein the step 
of replacing high intensity regions further comprises 
the step of replacing each one of said respective 

is sub-window having a quantity of high intensity pix- 
els greater than said predetermined intensity per- 
centage with the respective nearest sub-window 
having a quantity of high intensity pixels less than 
said preselected intensity percentage. 

20 

14. The method as recited in claim 1 3, wherein the step 
of replacing high intensity regions further comprises 
the step of selecting said predetermined intensity 
percentage to be about fifty percent. 

25 

15. The method as recited in claim 14, wherein the step 
of converting said modified x-ray image to the fre- 
quency domain further comprises the step of re- 
moving ringing artifacts. 

30 

16. The method as recited in claim 15, wherein the step 
of removing ringing artifacts comprises the step of 
replacing each adjacent window to a respective 
modified window with an equivalent modified win- 

35 dow to generate a super-window. 

17. The method as recited in claim 16, wherein the step 
of converting said modified image to the frequency 
domain further comprises the step of converting 

40 said super-window to the Fourier transform domain 
so as to generate at least one object data spectral 
component and at least one grid line spectral com- 
ponent that corresponds to said modified x-ray im- 
age. 

45 

18. The method as recited in claim 1 7, wherein the step 
of removing grid line artifacts comprise the steps of: 

selecting a respective one of said grid line arti- 
so fact spectral components being situated adja- 

cent to said at least one object data spectral 
component that is substantially lower in ampli- 
tude than said at least one grid lino artifact com- 
ponent; and 

55 

replacing said at least one grid line artifact com- 
ponent with a replacement spectral component 
having an amplitude substantially equal to the 


45 

18. 


so 
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average magnitude of said at least one obj ct 
data spectral component. 

19. The method as recited in claim 18, where the step 
of removing grid line artifacts further comprises the 
step of selecting said grid tine spectral component 
having a respective magnitude about three stand- 
ard deviations greater than the average magnitude 
of said at least one object data spectral component. 

20. The method as recited in claim 1 9, wherein the step 
of restoring said modified image comprises the step 
of discarding each of said equivalent modified win- 
dows from said super-window. 

21 . The method as recited in claim 20, wherein the step 
of restoring said modified image further comprises 
the step of restoring each of said edgy regions 
which were replaced. 

22. The method as recited in claim 21 , wherein the step 
of restoring said modified image further comprises 
the step of restoring each of said high intensity re- 
gions which were replaced. 

23. The method as recited in claim 1 , wherein the step 
of dividing said x-ray image into windows further 
comprises the step of calculating a pixel mean value 
of each respective window in said x-ray image. 


25. The method as recited in claim 1 , wherein the step 
of dividing said x-ray image into a predetermined 
number of windows comprises the step of dividing 
said x-ray image into sixty-tour windows. 


10 


1S 


20 


2S 


30 


24. The method as recited in claim 23. wherein the step 
of restoring said modified image further comprises 
the step of equalizing the pixel mean value of each 
said respective window by adding the difference be- 
tween the pixel mean value of each said respective 35 
restored modified window and the pixel mean value 
of each respective said window so as to equalize 
the mean value of said respective restored window 
to the mean value of said respective window. 


40 


45 
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